hep-ph/0303192 



Global Analysis of a Supersymmetric Pati-Salam Model 

T. Blazek*, S. F. King and J. K. Parry 

Department of Physics and Astronomy, University of Southampton 
Southampton, SO 17 IB J, U.K 
Abstract 

We perform a complete global phenomenological analysis of a realistic string-inspired 
model based on the supersymmetric Pati-Salam SU(4) x SU{2)l x SU(2)r gauge 
group supplemented by a U(l) family symmetry, and present predictions for all ob- 
servables including muon g — 2, r — > fij, and the CHOOZ angle. Our analysis 
demonstrates the compatibility of such a model with all laboratory data including 
charged fermion masses and mixing angles, LMA MSW and atmospheric neutrino 
masses and mixing angles, and b —>■ 37, allowing for small deviations from third fam- 
ily Yukawa unification. We show that in such models the squark and slepton masses 
may be rather light compared to similar models with exact Yukawa unification. 
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1 Introduction 



The origin of fermion masses and mixing angles represents a challenge faced by theo- 
rists for a long time. In the post-SuperKamiokande era this puzzle has become more 
intriguing than ever before. SuperKamiokande evidence for atmospheric neutrino os- 
cillations [1] has taught us that neutrino masses are non-zero and furthermore that 
the 23 mixing angle is almost maximal. More recently SNO [2] and KamLAND [3] 
experimemts have confirmed the matter enhanced Large Mixing Angle (LM A) solu- 
tion to the solar neutrino problem [4]. In this work we assume that the smallness 
of neutrino masses can be explained by the see-saw mechanism involving very heavy 
right-handed neutrino states, and that the see-saw mechanism is implemented using 
single right-handed neutrino dominance [5] which can explain in a natural way the 
coexistence of large neutrino mixing angles with a mass hierarchy. It then becomes a 
flavour problem to fit together the neutrino mass puzzle with the pieces provided by 
the long-known pattern of quark and charged lepton masses. 

The flavour problem cannot be fully addressed without unification. However, 
unification has its own challenges. These include the unification of gauge couplings 
and third family Yukawa couplings and the introduction of super symmetry. It is 
well known that supersymmetry facilitates gauge coupling unification, stabilises the 
hierarchy between the high energy scale and the weak scale, and allows a radiative 
mechanism of electroweak symmetry breaking. Within the natural framework of su- 
persymmetric unification, the larger high energy gauge group in turn increases the 
predictive power of the theory in the flavour sector, for example by leading to group 
theoretical mass relations between quark and lepton masses of the same family. Re- 
lations between quarks and leptons of different families require an additional family 
symmetry, however. In this way it becomes possible to address both the flavour 
problem and the unification problem, within a single framework. Having defined 
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the framework, it is by no means guaranteed that models exist which satisfy all 
the phenomenological constraints provided by current data, and comply with all the 
theoretical requirements such as successful electroweak symmetry breaking, and ap- 
proximate gauge and Yukawa unification, while reproducing the known observables. 
It is therefore important to know that at least some models exist which satisfy all 
the constraints, as an existence proof that such a proceedure can be implemented 
consistently. 

In this paper we shall study a particular example of a complete supersymmetric 
unified model of flavour, based on the Pati-Salam SU(A) x SU(2) L x SU(2) R gauge 
group [6] extended by an additional U(l) family symmetry. Accepting minimality as 
a model building principle this group has the following nice features: it establishes 
the third family Yukawa unification, places the right-handed neutrinos into non-trivial 
multiplets and does not introduce unwanted exotic states in the multiplets containing 
the Standard Model fermions and two Higgs doublets required by its SUSY extension. 
The Pati-Salam group can emerge from a simple gauge group like SO (10) or E(6). 
However, from a string theory perspective, it is not necessary in order to achieve 
unification that there should be a unified field theory based on a simple group. A 
partially unified gauge group can equally well emerge directly from string theory, 
and in the case of the Pati-Salam gauge group this possibility has been explored 
extensively both in the case of weakly coupled fermionic string theories [7] and in the 
case of type I strings with D-branes [8] . 

Although models based on the Pati-Salam gauge group have been extensively ex- 
amined, there is currently no complete up to date phenomenological study of this 
model in the literature. For instance [9] investigated constraints from Lepton Flavour 
Violation(LFV) in a Pati-Salam model with small neutrino mixing angles. Subse- 
quently a Pati-Salam model was proposed [10], using single right-handed neutrino 
dominance [5] to achieve naturally large neutrino mixing angles, but the question of 
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LFV was not readdressed and it was later shown [11] that the // — > branching 
ratio is too large. Moreover, only the negative /i parameter was considered in [10, 12] 
which is currently disfavoured by the muon g — 2. In other works such as [13] and 
[14] the neutrino sector is absent all together. The complete lepton sector is studied 
in great detail in a global analysis in [11], but the quark mass matrices used [10] 
were obtained for the opposite sign of ji, and the analysis gives imperfect fits for the 
branching ratio BR{b — > 57) or b quark mass mj, which both get potentially signifi- 
cant contributions from SUSY loops proportional to fi. To summarise, a completely 
phenomenologically acceptable supersymmetric Pati-Salam model does not currently 
exist in the literature. This illustrates the broader point that while many models exist 
in the literature, it is less common for the analysis of any such model to be complete. 

In this paper, then, we shall construct a "4221" model, following the approach 
of [10], and demonstrate its phenomenological viability. The model has approximate 
third family Yukawa unification perturbed by higher order terms and assumes non- 
universal soft Higgs masses. To demonstrate the viability of such a model, we perform 
a top down global analysis of the parameter space carried out on 24 observables. In 
the leptonic sector the observables include the muon g — 2 and solar and atmospheric 
neutrino data. A complete list of observables and their a values, which are used to 
calculate the x 2 function can be found in Table 2. In the analysis we ensure that the 
upper limits on the branching ratio for the lepton flavour violating processes r — > /ry, 
fi — > e7 and r — > &y are not exceeded as well as the limit on the 13 neutrino mixing 
angle. In addition to this an experimental lower bound on each sparticle mass was 
imposed. In particular, the most constraining are: the LEP limits on the charged 
SUSY masses (m^,m f > 105GeV), the CDF limit on the mass of the CP odd Higgs 
state (m A o > 105-110 GeV, valid for tan (3 50) [15], and the requirement that the 
lightest SUSY particle should be neutral. 

Having constructed the model and demonstrated its phenomenological viability, 
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we then discuss the following three aspects of the model in more detail: 



• The first such aspect, as first pointed out in [11], is lepton flavour violation 
arising from the large 23 neutrino mixing through a neutrino Yukawa texture 
of the form 



Due to large tan (3 additional features emerge when studying correlations among 
observables like BR(6 — * sj), BR(r — > /ry) and muon g — 2. Most notably, 
two distinct minima are found with similar x 2 values for the best fits. These 
conclusions are new since study [11] did not investigate a complete model and 
all other previous works did not involve global analysis. 

• The preference for positive //, given by the sign of the muon g — 2 discrepancy, 
implies positive gluino corrections to thus leading to difficulties in obtaining 
t — b — t Yukawa unification. Hence a second focus of the present work is to 
study the required deviation from third family Yukawa unification in the best 



• Thirdly we focus on the effects of future experimental advances, in the form of 
direct Higgs searches, a lepton flavour violating r — > /ry measurement and a 
refinement of the muon g — 2 discepancy, upon our global fits, indicating how 
further experimental progress in these areas will constrain the parameter space 
of the model. 

The remainder of the paper is layed out as follows. Section 2 briefly reviews 

our construction of a string-inspired Pati-Salam model. Section 3 contains a brief 

description of the numerical technique used in the analysis. A discussion of our main 

results can be found in section 4, with concluding remarks in section 5. 
lr This was also recently studied from a somewhat different point of view in [13]. 
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2 The Model 



The model considered in this paper is based on the Pati-Salam gauge group [6], 
supplemented by a U(l) family symmetry, 

SU(4) x SU(2) L x SV(2) R x U{1) (2) 

The left-handed quarks and leptons are accommodated in the following representa- 
tions, 

f- = ( 4,2,i) = (£ :l % ;)' 0) 

— ■ - - / d R d B <F e+ V 

(4,1,2)= \ % % \ (4) 



where a = 1 ... 4 is an SU(4) index, a,x = 1, 2 are SU(2) LjR indices, and % — 1 . . . 3 
is a family index. The Higgs fields are contained in the following representations, 

(where /ii and h 2 are the low energy Higgs superfields associated with the MSSM.) 
The two heavy Higgs representations are [7] 

!T* = (4,l,2)=(>f "| "I ^ ) (6) 
V e # / 



and 



ff- = (J,l,5)=| | | | £ )• (7) 



The Higgs fields are assumed to develop VEVs, 

< H >=< u H >~ M GUT , < H >=< u H >~ M GUT (8) 
leading to the symmetry breaking at Mqut 

SU(4) <g> SU(2) L <g> SU(2) R — ► SU(3) C <g> SU(2) L <g> U(l) y (9) 
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in the usual notation. Under the symmetry breaking in Eq.9, the Higgs field h in 
Eq.5 splits into two Higgs doublets hi, hi whose neutral components subsequently 
develop weak scale VEVs, 

< h\ >= vi, < h G 2 >= v 2 (10) 

with tan/3 = v 2 /vi. 

To construct the quark and lepton mass matrices we make use of non-renormalisable 
operators [16] of the form: 

») im(m(m m (±Y". m 



\M 2 J \M 2 J \M / 

The 9 fields are Pati-Salam singlets which carry U(l) family charge and develop VEVs 
which break the U(l) family symmetry. They are required to be present in the op- 
erators above to balance the charge of the invariant operators. After the H and 9 
fields acquire VEVs, they generate a hierarchy in i) effective Yukawa couplings and 
ii) Majorana masses. These operators are assumed to originate from additional inter- 
actions at the scale M > M GUT . The value of the powers p^ and are determined 
by the assignment of U(l) charges, with Xg = —1 then = (X F i + X-pj + X h ) and 
q l3 = (X T +X yJ +X h ). 

The contribution to the third family Yukawa coupling is assumed to be only from 
the renormalisable operator with n = p = leading to Yukawa unification. The 
contribution of an operator, with a given power n, to the matrices Yf =u ^, u ,e, Mrr 
is determined by the relevant Clebsch factors coming from the gauge contractions 
within that operator. A list of Clebsch factors for all n — 1 operators can be found 
in the appendix of [10]. These Clebsch factors give zeros for some matrices and not 
for others, hence a choice of operators can be made such that a large 23 entry can be 
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given to Y v and not Y u ,d,e- We shall write, 

<HxH> <9> 
6 = Af2 = °- 22 > e = l^ = - 22 ' ( 13 ) 

then we can identify S with mass splitting within generations and e with splitting 

between generations. 

Our choice of U(l) charges are as in [10], and can be summarised as Xpi = (1, 0, 0), 
Xpi = (4, 2, 0), X h — 0, X H = and Xjj = 0. This fixes the powers of e in each 
entry of the Yukawa matrix, but does not specify the complete operator. The Yukawa 

couplings are specified by the choice of operators, 

/ (a u O R + < 1 0" y )e 5 (a 12 T + a' 12 0' F )e 3 (a' 13 0' c )e 
Y f (M GUT )=\ (a 21 G y (a 22 O w + a' 22 0' s )e 2 (a^O 1 + a' 23 0' w ) 

V (aaiO^e 4 (a 32 M + a' 32 0' K )e 2 a 33 

(14) 

The operators are defined in [10], although the selection of operators here is different 
from that paper. The notation is such that O, O' and O" are n — 1, n — 2 and 
(highly small) n = 3 operators respectively where n refers to the powers of (HH), 
thus 2 

O ~ (M) ~ 5, C ~ (HH) 2 ~ 5 2 , 0" ~ (HH) 3 ~ <5 3 . (15) 

The order unity coefficients Oy,a^,a^ multiply the operators 0,0',0" in the r/' 
position. The Majorana operators are assumed to arise from an m = operator 
in the 33 position and m — 1 operators elsewhere. The operator choice in Eq.14 
leads to the quark and lepton mass matrices in Table 1. For example the Clebsch 
coeficients from the leading O w operator in the 22 positon gives the ratio : 1 : 3 in 
the Y\j,d,e matrices. This ratio along with subleading corrections provides the correct 
m c : m s : m M ratio. 

In the neutrino sector the matrices in Table 1 satisfy the condition of sequential 

dominance [5] in which a neutrino mass hierarchy naturally results with the dominant 

2 Thc n = 3 operators can, to a very good approximation, be neglected. Their inclusion here 
serves only to fill the 11 entries of the Y u v Yukawa matrices, thereby ensuring (for example) that 
the up quark is given a very small mass. 



Y U (M GUT ) = 



( V^a'l^e 5 V2a[ 2 5 2 e 3 



V5 



'13 



5 2 e\ 



343 
670 



a' 22 S 2 e 2 



a' 32 S 2 e 2 



r t a 3 3 ) 



Y d (M GUT ) 



Y e (M ( 



GUT 



Y V {M ( 



GUT 



M RR (M ( 



GUT 



M 
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( - an5e 5 
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—j= a 2 i5e 4 
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V 
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A 12 5e 6 
A 13 5e* 



-V2a'S 2 e 3 



a 22 5e s 



1372 
"670" 



a' 22 S 2 e 2 
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-3\l-a 22 5e 2 + \/—a 22 5 2 € 2 



772 



670 



V2a 32 Se 2 



2 ai 2 Se 3 
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193 
670 

6 



a' 22 5 2 e 2 2 a 23 S 



a' 32 5 2 e 2 r v a 33 



A 12 5e e A 13 5e 4 \ 
A 22 5e A A 23 5e 2 



A 23 5e 2 A 



/ 



a' 23 5 2 



r b a 33 J 



-3\l^a> 23 5 2 



a 33 



Table 1: The quark and lepton Yukawa matrices and neutrino Majorana mass ma- 
trix as used in the analysis. In our numerical analysis we set Mr = 3 • 10 14 GeV. 
The Yukawa matrices follow from Eq.14 and the Clebsch factors arising from each 
operator are shown numerically above. Clebsch zeroes play an important part in 
suppressing the leading operator contribution in a particular element of the matrix, 
or in simply giving a zero if all the operators are suppressed. The Clebsch coefficients 
in the Majorana sector are set equal to unity with Aij being independent order unity 
coefficients. 
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third right-handed neutrino being mainly responsible for the atmospheric neutrino 
mass, and the sub-dominant second right-handed neutrino being mainly responsible 
for the solar neutrino mass. The atmospheric mixing angle is then determined approx- 
imately as a ratio of Y 23 : Y 33 , and the solar mixing angle is determined by a ratio of 
Yj} 2 to a linear combination of Y 22 and Y 32 , while the CHOOZ angle is determined by 
a more complicated formula [17]. Note that the dominant right-handed neutrino in 
this model is the heaviest one, corresponding to heavy sequential dominance (HSD) 
and LFV has been considered in general in this class of models [18]. 

In the previous analysis [11] the matrix elements, F e 12 , Y^ 3 were suppressed ar- 
tificially to keep BR(fi — ► ej) within its experimental limit without substantially 
changing the predictions of fermion masses and mixings. In this new analysis we 
have built this suppression into the model with our new choice of operators, whose 
Clebsch coefficients give zeros in the desired matrix elements as can be seen in Table 1. 
This can be understood analytically from [18]. 

The subleading operators in the 33 position are not shown explicitly, but are 
expected to lead to significant deviations from exact Yukawa unification. This effect 
is parametrised by the ratios 

_ Y U (M GUT ) 33 _ Ycl(Mgut) 33 _ Y U (M GUT ) 33 
n ~ Y e (M GUT )^ n ~ Y e (M GUT ) 33 > ru ~ Y e (M GUT ) 33 - [ } 

3 Numerical Procedure 

In our numerical analysis we have adopted a complete top-down approach [19]. At 
the GUT scale the MSSM gauge couplings are related to the GUT scale couplings 
as a 2 L = a i = ®gut and a 3 = a G uT(l + £3), where e 3 sums up the effects of GUT 
scale threshold corrections. The particular choice of the Yukawa couplings, Table 1, 
follows from the higher dimensional operators in Eq. (14) as the latter are matched 
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to the MSSM lagrangian. The parameters 

M GUT , olcut-i e 3 , 6, e, a's and A's, r t , r 6 , rv, 

(17) 

M1/2, ^4o, A*, -B/i, ra|, m^, ra£ and £> 2 
are then defined by the boundary conditions at the GUT scale. They parametrise 
the imprint of a complete Pati-Salam theory together with the SUSY sector (second 
line) on the MSSM and stand for the input in case this theory is not fully known. In 
the SUSY sector, the soft SUSY breaking parameters are for simplicity introduced at 
the same scale. The gaugino masses are assumed universal (equal to M1/2) and so do 
the trilinear couplings: A4, = AqY^, for % = u, d, e, v. The soft scalar masses of the 
MSSM superfields include the D terms from the breaking of the Pati-Salam gauge 
group [12] 

2 2 1 2 r^2 

m Q = m F + g 4 D 

m ln = ra\-{g\-2gl R )D 2 

m X = m^-(gl + 2g 2 2R )D 2 

m 2 L = m 2 F -3glD 2 (18) 

m l R = rn^+(3g 2 4 -2g 2 2R )D 2 

m l R = m^+(3g 2 4 + 2g 2 2R )D 2 

m H u = rn 2 h -2g 2 2R D 2 

™ 2 Hd = ™l + Zg 2 2R D 2 . 

As D 2 = I ( \H U \ 2 — \H U \ 2 ) [12] it is possible for this quantity to be both positive 
and negative. 

We now describe minor simplifications to the input in (17) which were assumed 
in the actual numerical analysis. We have kept equality between the two order pa- 
rameters 5 and e as in Eq.(13) and the soft SUSY breaking scalar masses m F and 
m-p have been held equal to each other as well. Furthermore we exploited the fact 
that determining h(Mgut) and Bfi at the GUT scale is equivalent to determining the 
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low energy values fi(M z ) and tan/3, respectively. Thus instead of (17) our numerical 
analysis uses 

M GUT , a GUT , e 3 , 5, a's and A's, r t , r b , r„, 

(19) 

Mu2, A , fi(Mz), tan/3, m 2 F , m\ and D 2 
as input parameters. The top down approach implies that we can freely vary or 

hold fixed any one of them and then investigate the fit properties. This is one of 

the advantages of doing the analysis top down. For example, in more traditional 

bottom up approaches it is difficult to control the size of the dimensionless GUT 

scale parameters. One usually sets up a sample of randomly scattered points and 

then searches through it to identify a sub-sample with physically interesting GUT 

scale properties. In our case we can set up the interesting GUT relations explicitly 

right at the start — as we have done for instance in section 4.3 where the fits with r& 

and r t approaching unity are studied. 

We note that taking advantage of the top-down approach we kept S = 0.22, r v = 1, 
A = and tan j3 = 50 fixed throughout the analysis. This effectively increases the 
degrees of freedom of the global analysis and can provide a reference point if further 
analysis is required in the future. We also kept the \i parameter at scale M z fixed to 
two different values as is explained below. 

Two-loop RGEs for the dimensionless couplings and one-loop RGEs for the di- 
mensionful couplings were used to run all couplings down to the scale M 3R where the 
heaviest right-handed neutrino decoupled from the RGEs. Similar steps were taken 
for the lighter M 2 r and M 1R scales, and finally with all three right-handed neutrinos 
decoupled the solutions for the MSSM couplings and spectra were computed at the 
Z scale. This includes full one loop SUSY threshold corrections to the fermion mass 
matrices and all Higgs masses while the sparticle masses are obtained at tree level. 

rrih and D in Eqs. (19) were varied to optimise radiative electroweak symmetry 
breaking (REWSB), which was checked at one loop with the leading m\ and 
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corrections included following the effective potential method in [20]. We note that 
as tan/3 determines the Higgs bilinear parameter B/i, there is a redundancy in our 
procedure since two input parameters, rrih and D, determine one condition for the 
Higgs VEV of 246 GeV. This approach enabled us to control the /j, parameter and we 
explored regions with \i low (// = f20GeV) and high (// = 300GeV) 3 . 

An experimental lower bound on each sparticle mass was imposed. In particular, 
the most constraining are: the LEP limits on the charged SUSY masses (m^± , > 
105GeV), the CDF limit on the mass of the CP odd Higgs state {m A o > 105-110 GeV, 
valid for tan/3 w 50) [15], and the requirement that the lightest SUSY particle should 
be neutral. Finally, the \ 2 function ^2(X* h — X^ xp ) 2 /af) is evaluated based on 
the agreement between the theoretical predictions and 24 experimental observables 
collected in Table 2. In addition to the constraints listed above and in [11], we make 
a full analysis of the quark sector mass and mixings, in particular we have included 
the important constraint set by BR(6 — > sj). 



3 For tan/3 as large as 50, fi » 300GeV leads to too large SUSY threshold corrections to the 
masses of the third generation fermions r and b unless the sparticles in the loop have masses well 
above the 1 TcV region. [21, 19] 
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Observable Mean Oi 





1/137.036 


7.30 • 10~ 6 


G ii 


1.16639 • 10~ 5 


1.12 • 10~ 7 




0.1181 


0.0020 


M t 


174.3 


5.1 


m h {m h ) 


4.20 


0.20 


M b - M c 


3.4 


0.2 


m s (2GeV) 


0.110 


0.035 


9 9 \ / 9 
m d ~ m u)/ m s 


2.03 • 10~ 3 


2.0 • 10~ 4 


m d /m s 


0.05 


0.015 


M T 


1.777 


1.8 • 10~ 3 


M„ 


0.106 


1.1 • 10~ 4 


- LV± e 


5.11 • 10~ 4 


5.1 • 10~ 7 


| v us \ 


0.2196 


0023 


\V,h\ 

\ v CO \ 


0.0402 


0.003 


\v ub \/\v cb \ 


0.09 


0.02 


M z 


91.1882 


0.091 


M w 


80.419 


0.08 


Pnew 


-0.0002 


0.0011 


BR{b -> si) 


3.47- 10~ 4 


0.45 • 10~ 4 


£<V NEW 


34.7.6 • 10~ 10 


11 - 10 10 


Am ATM 


2.5 • 10~ 3 


0.8- 10~ 3 


sin 29 1 atm 


0.99 


0.06 


Am2 SOL 


7.0- 10~ 5 


3 • 10~ 5 


sin 2 26 sol 


0.8 


0.09 



Table 2: Table of observables and a values which are used to calculate the \ 2 that enables 
best fit regions to be determined via minimisation. 
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4 Results and Discussion 



The numerical results from the global analysis are presented in the form of contour 
plots in the (mp, M1/2) plane and are produced for two different values of the mu 
parameter \x = 120 GeV and /x = 300 GeV. Before we address the details we would 
like to discuss two different viewpoints of our analysis, namely the flavour sector on 
the one hand and the unification sector of the other hand. In our discussion we would 
like to distinguish between the two viewpoints. The main distinction is that in the 
MSSM analysis the flavour parameters (with the exception of a 33 ) and Aij can 
be considered fixed at unity or at a value of order unity. Up to a 33 which enters the 
large Yukawa couplings their exact values do not affect the fit of the SUSY spectra 
or SUSY-related observables like the muon g — 2 or branching ratio b — > 57. They 
neither perturb gauge coupling unification nor change the running of the large Yukawa 
couplings. This means that the discussion of our results is naturally split into a part 
dealing with the flavour structure of the Pati-Salam model where the variation of the 
coefficients of the higher dimensional operators matters, and a part where the MSSM 
analysis is presented and the conclusions do not depend on the variation of the a and 
A parameters (up to a 33 ). 

Concerning the flavour sector, our results can be used to show how well the model, 
i.e. the set of higher dimensional operators specified by Eq. (14), describes the ob- 
served fermion masses and mixings. Taking this viewpoint all parameters listed in 
(19) represent the input of the analysis. The results in either of the four panels in Fig- 
ure 1 show that the model gives a very good agreement with the data. The minimum 
of the total x 2 is less than unity obtained for \i = 120 GeV in the upper left panel. 
This means that it is possible to fit every observable to better than a 1 er accuracy. 

Concerning the unification sector, the conclusions are much stronger as much fewer 
number of the input parameters enters effectively after the a's and A's decouple from 
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the analysis. Indeed, the set of the effective input parameters in this sector is reduced 
to 

M GUT , a GUT , e 3 , a 33 , r t , r b , 

(20) 

M1/2, rrip, ml and D 2 . 
With this input the low energy Higgs and SUSY spectra are determined. The con- 
ventional present-day observables include a em, G>, a s {Mz), M t , rrib(mb), M T , Mz, 
Mw, Pnew, BR(b — > sj) and 5a M . The a^- and ^4^- input parameters are all of order 
one and their exact values are always adjusted to fit the first two generation masses 
and mixings well while these variations do not affect the fit of the observables listed 
above. 

We study many details of the MSSM analysis, in particular the dependence of the 
fit on mp and My 2 , best fit results for muon g — 2, and BR{r — > /ry) predictions. 
The numerical results also contain studies of a deviation from Yukawa unification and 
a future measurement of BR(t — > fxj). The effect of a change to the present muon 
g — 2 discrepancy was studied and also the effect of future direct Higgs searches, the 
results of which can also be found at the end of the paper. 

From our global analysis we found that there are two x 2 minima as shown in 
Figure 1. In this model there are two conditions and three free variables, m\, D 2 , 
and Bfi, for electroweak symmetry breaking to be achieved. The two minima hence 
are independent solutions to these conditions. Minimum A has D 2 negative and 
smaller rrih, Bp. Minimum B on the other hand has D 2 positive and larger m^, Bp. 
The relative size of B\i results in a different Higgs spectrum, particularly the CP odd 
pseudoscalar Higgs, A , which will be lighter for minimum A and heavier for minimum 
B. The difference between the sign of D 2 , which contributes to the soft scalar masses 
as shown in Eq. (26), means that minimum B will have lighter right squarks and left 
sleptons, along with heavier left squarks and right sleptons, than minimum A. This 
difference in sign of D 2 has some interesting phenomenological consequences for the 
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U.oU 


^21 


1.0/ 


1.0 / 


1.0 1 


1. 1 


^22 


1 1 3 




1 1 3 


1 1 3 


u 23 


0.98 


0.89 


1.05 


0.85 


031 


-0.20 


-0.21 


-0.20 


-0.28 


032 


2.18 


2.08 


2.32 


2.53 


«12 


0.77 


0.77 


0.71 


0.71 


«i 3 


0.60 


0.53 


0.46 


0.46 


°22 


0.66 


0.66 


0.64 


0.62 


°23 


0.41 


0.40 


0.36 


0.36 


a 32 


1.16 


1.80 


1.56 


1.72 


<1 


0.32 


0.278 


0.20 


0.23 


^11 


0.63 


0.94 


0.63 


0.94 


^12 


0.74 


0.48 


0.69 


0.52 


^13 


1.75 


2.10 


1.73 


2.04 


^22 


0.97 


0.52 


0.93 


0.55 


^23 


2.49 


1.79 


2.23 


1.91 


^33 


1.97 


1.88 


1.97 


1.88 



Table 3: Tables of inputs for the best fit points for each of the global \ 2 minima with 
H = 120 and 300 GeV. 
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Outputs 




/x = 120 GeV 


/x = 300 GeV 




Min A 


Min B 


Min A 


Min B 


m,A° 


102 


818 


102 


822 


m h o 


106 


114 


106 


114 


m H o 


112 


891 


113 


888 


m H + 


136 


861 


135 


861 


Mi 


186 


270 


186 


271 


M 2 


371 


537 


371 


537 


M 3 


1175 


1671 


1175 


1671 


M xt 


114 


117 


272 


290 


M xi 
M Nl 


390 


549 


408 


554 


98 


107 


179 


249 


M N2 


130 


127 


277 


305 


M N 3 


198 


278 


307 


311 




390 


549 


408 


554 


M Q 12 


1166 


1679 


1159 


1673 


Mq 


979 


1345 


960 


1356 




1131 


1623 


1124 


1617 


M u 3 


798 


1147 


805 


1160 


M D 12 

Mr* 


1182 
923 


1510 
1192 


1204 
1044 


1529 
1251 


M Ll 


673 


611 


715 


656 


M~ L2 


665 


595 


707 


644 


M U 


580 


334 


638 


425 


M El 


496 


766 


473 


752 


M E2 


495 


765 


473 


751 


Me, 


201 


370 


188 


325 


t -> /x 7 


2 • 10" 7 


3 • 10" 6 


8 • 10" 8 


5 • 10" 7 


r — > e7 


1 • 10~ 14 


3 • 10~ 13 


6 • HT 15 


5- 10~ 14 


// -> e7 


3 • 10- 14 


1 • 10" 13 


1 • 10" 14 


3 • 10" 14 


sin 6*i3 


0.053 


0.078 


0.037 


0.10 


sin(/3 — a) 


0.22 


1.0 


0.15 


1.0 


cos(/9 — a) 


-0.98 


0.0 


-0.99 


0.0 



Table 4: Tables of outputs for the best fit points for each of the global \ 2 minima with 
H = 120 and 300 GeV. The input parameters are as defined in Table 3. 
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two minima which will now be discussed. 

The upper and lower plots shown in Figure 1, display the \ 2 contours for these 
two minima. Each of the figures display results for both /i — 120 and 300 GeV in 
the relative left and right positions. The contours in Figure 1 are bounded from the 
lower Trip region due to the lightest stau becoming the LSP and from the lower Mi/ 2 
region due to an increasing x 2 penalty coming from BR{b — > 57). 

The upper minima of Figure 1, Minima A, have a preferred region in the lower 
(m F , M1/2) plane, with M 1/2 = 400 - 500 GeV and m F = 500 - 700 GeV. The 
lower minima of Figure 1, Minima B, have their preferred region nearer Mi/ 2 = 
550 — 650 GeV and mp = 600 — 800 GeV. A list of inputs and outputs for the best 
fit point in each minimum can be found in Tables 3 and 4. The Higgs masses and 
CP even Higgs mixings found for minimum A in Table 4 are discussed in detail in 
section 4.4. 

4.1 Muon g - 2 

Figure 2, shows contour plots for the SUSY contributions towards the muon g — 2. 
Both minimum A and minimum B (upper and lower plots respectively) give good fits 
to the present discrepancy between experiment and Standard Model prediction. As 
expected, a larger contribution to the muon (7 — 2 is obtained in the lower left corner 
of the (m F ,M 1 / 2 ) plane where the SUSY spectrum is lightest and decreases as we 
move towards a heavier spectrum in the top right corner. It is also clear that for any 
one point in the (mp, Mi/ 2 ) plane, minimum B gives a larger contribution than the 
corresponding point in minimum A. This relative enhancement can be ascribed to the 
dominant chargino-sneutrino diagram via the presence of a lighter muon sneutrino 
for the case of minimum B, as can be seen in Figure 3. 

The present muon g — 2 discrepancy lies at 34 x 10~ 10 but over the past 12 months 
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it has varied from a 1.5a to 3a effect. Also the size of the present discrepancy depends 
on the experimental data used in the calculation of the Standard model prediction. 
The value we have used throughout our analysis [23] makes use of e + e~ data. On 
the other hand it is possible to do the same calculation making use if r decay data 
[24], which gives a lower discrepancy of 9.4 x 10~ 10 . As a result we think it worth 
while looking into how our best fit regions would change if a lower discrepancy was 
assumed. For simplicity we took 3 points in the (mp, My 2 ) plane of minimum A with 
ji = 120 GeV and gradually changed the g — 2 discrepancy from 34 x 10~ 10 down 
to 0. The results are presented in Figure 4 as a plot of x 2 against the muon g — 2 
discrepancy, a^ ew . 

With the discrepancy held at 34 x 1CT 10 the best fit point is near Mi/ 2 = 450 GeV 
and m F = 550 GeV. Following the curve corresponding to this point in parameter 
space, we can see that as the muon g — 2 discrepancy is lowered the x 2 gradually 
increased. Therefore the best fit point has moved in the positive M 1 / 2 , ™>f direction. 
Looking at the two further curves in Figure 4 we can see that if a^ ew ~ 16 x 10~ 10 
then the best fit point would move nearer My 2 = 550 GeV and m F = 650 GeV. One 
particular point of interest is a^ ew = 9.4 x 10~ 10 , the value for the discrepancy as given 
by the Standard Model prediction from r decay data. If we make an approximation, 
based on the curves in Figure 4, we can say that the best fit point, for a^ ew = 
9.4 x HT 10 , would be in the region M 1/2 = 550 - 700 GeV, m F = 650 - 700 GeV. 

4.2 r — > ^7 

Figure 7 displays contours for the quantity BR(r — > 717) for both minima with /i = 120 
and 300 GeV as labelled. The general pattern of the contours show larger branching 
ratio for lighter SUSY spectrum and smaller branching ratio for heavier spectrum. 
This pattern is not strictly obeyed in the bottom left panel which shows results for 
minimum B with /i = 120 GeV. The reason for this is that our numerical procedure 
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adds a large penalty x 2 contribution for a r — > /ry branching ratio larger than the 
BaBar limit of 2.0 x 10~ 6 . Looking at the bottom left panel in Figure 7 we would 
expect the branching ratio to exceed the present limit as we go to a lighter spectrum. 
The result of adding this penalty x 2 is to numerically force an alternative solution to 
be found which gives lower branching ratio and disrupts the pattern. Recalculation 
of this region of parameter space without the additional x 2 penalty does indeed yield 
values of BR(r — > /ry) as large as 6.0 x 10~ 6 , these points would therefore follow the 
expected contour pattern but are clearly experimentally excluded. 

Looking at Figure 7, the branching ratio for minimum A with \i = 120 and 
300 GeV is well below the present experimental bound. On the other hand, the 
branching ratio for minimum B, with /i = 120 GeV Figure 7, is right at the present 
90% confidence level bound of 2.0 x 10~ 6 [25]. For fi = 300 GeV minimum B gives a 
branching ratio in the range, 0.1 — 0.2 x 10~ 6 , just below the present bound. With 
BaBar expected to search as far as BR(r — > /ry) ~ 10" 8 over the next 5 years this 
certainly provides a means of distinguishing the two minima. 

4.3 Deviations from Yukawa Unification 

The plots shown in Figure 8 show contour lines for r b = Y b /Y T and those in Figure 9, 
show contour lines for r t = Y t /Y T in the best fits over the (mj? , M 1( / 2 ) plane. These 
parameters allow the deviation of the top, bottom and tau Yukawa couplings away 
from unification(r fe = r t = 1). Both parameters show significant dependence upon m F 
and weak dependence upon Mi/ 2 , with increasing r t) b values as we move towards larger 
Trip. The plots show that the level of deviation from Yukawa unification required for a 
good x 2 St to be obtained is of the order of 20-35% in r b and 0-10% in r t . It is possible 
to account for this level of deviation through the presence of subleading operators, of 
the type mentioned in Eq. 15, in the 33 element of the Yukawa matrices. Hence the 
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33 element in Eq. 14 should read, 

Y 33 = a 33 + O + O' + . . . (21) 

where the operators O and O' are responsible for generating r t:b ^ 1. The 23 block of 
the neutrino Yukawa matrix has already shown us that a contribution to the Yukawa 
matrices from a subleading operator can actually be comparable to those from a 
leading operator. This occurs in the 23 element of the neutrino Yukawa matrix, where 
there is a contribution from the operator O 1 and the neutrino matrix is the only one 
that receives a non-zero Clebsch as can be seen in Table 1. This leads to the relative 
sizes of the elements Y V 2z ~ 0.44 and Y v iz ~ 1- A similar subleading contribution to 
the 33 element of the up and down quark Yukawa matrices could easily account for a 
deviation from third family Yukawa unification at the level discovered in our study. 

Here we do not study the region in the parameter space m F > 2 TeV, A « —2m F 
where the exact unification might work [26]. Instead, we carried out a study of the 
additional % 2 penalty incurred due to demanding exact Yukawa unification in the 
region tuf < 2 TeV and A = 0. Figure 10 shows the result as \ 2 contour plots in 
the r t — r b plane corresponding to the best fits. The three panels were obtained from 
three points in the (m^, Mi/ 2 ) plane and show that a very heavy penalty 5x 2 > 10 is 
paid when requiring exact Yukawa unification in this SUSY region. 

4.4 Future Higgs searches 

Figure 11 shows mass contours of the CP odd pseudoscalar Higgs, m^o. These plots 
show that for the Pseudoscalar Higgs mass minimum A prefers values approximately 
200 — 300 GeV lower than minimum B. In Figure 11 we see that for both fi = 120 and 
300 GeV, minimum A gives a very light pseudoscalar Higgs mass, m^o ~ 108 GeV, 
in the low M 1 / 2 , rap region. This is in fact the same region in which minimum A 
provides its lowest x 2 - m f &c t Table 4 shows that for the best fit point in minimum A 
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we have a pseudoscalar mass of 102 GeV and a light CP even mass of 106 GeV. With 
the Tevatron now taking data there is a high probability that the present lower bound 
on Higgs masses will be pushed higher. Hence we have undertaken a study of the 
effect this would have on our best fits. The plot in Figure 5 shows the increase in x 2 , 
for four points in the (mp, M1/2) plane of minimum A, due to an increase in the lower 
bound on the Higgs masses m^o, m h o. It clearly shows that all four of the points can 
accommodate an increase in the lower bound up to approximately 120 GeV, above 
this the x 2 increases sharply due the inability to accommodate such a large lower 
bound. 

The coupling of the light CP even Higgs, h°, to the Z boson is proportional to 
sin(/3 — a) and that of the heavy CP even Higgs, H°, is proportional to cos(/9 — a), 
where a is the mixing angle for the CP even Higgs states. In figure 6, which shows 
contours of sin(/3 — a) for points in minimum A, we see that in the low M1/2 region 
sin(/5 — a) is small and hence the Z couples dominantly to the heavier Higgs state H°, 
rather that the lighter h°. Therefore, in this region it is the heavier state, H°, which 
is the standard model like Higgs and so the LEP limit will apply to the larger m H o 
and not m h o . Table 4 shows that we have exactly this situation for the best fit points 
of minimum A where sin(/3 — a) ~ 0.2, therefore the standard model like Higgs is the 
heavier state H° for these points with a mass of 113 GeV. Assuming a 3 GeV error 
in our numerical calculation means we are compatible with the present LEP limit of 
114.4 GeV. 

We have checked that the light Higgs spectrum is consistent with the current limit 
on the rate for B s — > [22] 
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4.5 CHOOZ angle, 13 

Figure 12 shows scatter plots of sin 2 2#i 3 against Am\ tm for both minimum A and 
minimum B with \x = 120 and 300 GeV. Each point denotes values obtained from 
individual points in the (nip, Mi/ 2 ) plane, with points grouped according to the value 
of x 2 - These plots show that the Model can easily yield values of 13 that are within 
the present CHOOZ limit, #13 < 0.22. Each of the plots in Figure 12 shows that the 
best fit points, denoted by a + symbol, give a range of values of sin 2 2# 13 from, 10~ 6 
to 0.1. Although our results do not give any firm prediction for the value of #13, it 
can be seen that the model favours the region, 10~ 4 < sin# 13 < 0.1, just below the 
present CHOOZ limit. 

5 Summary and Conclusion 

We have performed a complete global phenomenological analysis of a realistic string- 
inspired model based on the supersymmetric Pati-Salam SU(A) x SU(2) L x SU(2) R 
gauge group supplemented by a U(l) family symmetry. Global contour plots in the 
(mp, My 2 ) plane have been presented in Figure 1, showing two x 2 minima. These 
two distinct minima differ numerically by the relative sign of the D-term. This gives 
interesting phenomenological differences between the two minima, notably one has 
BR(t — > (jlj) near the present limit and a heavy pseudoscalar Higgs m^o, while the 
other has BR(t — > /ry) well below the present bound but a light pseudoscalar Higgs 
m A o. Both minima give a good fit to the present muon g — 2 discrepancy over a large 
region of parameter space and give sin 2 2#i 3 over the range 10~ 5 — 0.1. Our best fit 
predictions for the superpartner masses for each of the two minima for two different 
fi values are summarised in Table 4. 

We emphasise again that our analysis really should be considered as consisting of 
two distinct parts, associated with flavour physics on the one hand and unification 
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and electroweak symmetry breaking on the other hand. For the flavour part, we 
have proposed a complete model in Table 1 which gives an accurate description of 
all fermion masses and mixing angles, including the LMA MSW neutrino solution. 
We have shown that improved limits on BR{r — > /ry) could begin to rule out one of 
our two minima. The conclusions on BR{t — > /ry) are applicable to a wide class of 
models which achieve approximate maximal atmospheric neutrino mixing via the see- 
saw mechanism in the MSSM with a large 23 entry in the neutrino Yukawa matrix. 
On the other hand BR(fi — > cy) is predicted to be about two orders of magnitude 
below the current limit, which is a consequence of the specific flavour structure of the 
model in Table 1. 

Regarding unification, the model predicts approximate third family Yukawa uni- 
fication and hence large tan/3 ~ 50. Electroweak symmetry breaking was achieved 
with the help of D-terms and non-universal soft Higgs mass, which allows small /i 
values. The property of exact Yukawa unification was relaxed throughout the analy- 
sis and it was found that a deviation of 20-35% for the bottom Yukawa coupling and 
0-10% for the top Yukawa coupling are required for a good fit to be obtained. We 
showed that relaxing Yukawa unification has the effect of allowing small values of the 
soft scalar mass mp, and lighter squark and slepton masses as a consequence. 

Further studies of the effects of future direct Higgs searches and a change to the 
present muon g-2 discrepancy are shown in Figures 4 and 5. We found that our best 
fit points, for the minima with lighter Higgs masses, can accommodate a lower bound 
on Higgs masses up to about 120 GeV. For these points the coupling of the lighter 
CP even Higgs state to the Z boson is suppressed, leaving the heavier of the two CP 
even states acting as the standard model like Higgs. 

In conclusion, we have constructed and analysed a complete supersymmetric Pati- 
Salam model which agrees with all laboratory observables and constraints. Using 
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a global analysis we identify the most preferred regions of the SUSY parameter 
space, and find a rather light superpartner spectrum corresponding to (mp, Mi/ 2 ) ~ 
(600, 600) (in GeV) well within reach of the LHC. 
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Figure 1: x 2 contour plot in the plane of (mF,Miy 2 )- The four plots, are obtained from 
the two minima, minimum A and minimum B with fi = 120 and 300 GeV as labelled. All 
points in the top left corner with approximately M ly / 2 > 700 GeV and mp > 700 GeV are 
unphysical due to the lightest stau becoming the LSP. 
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Figure 2: Muon g— 2 contour plot in the plane of (mp, M l7 / 2 )- The four plots, are obtained 
from the two minima, minimum A and minimum B with /j, = 120 and 300 GeV as labelled. 
All points to the left of the solid red line are unphysical due to the lightest stau becoming 
the LSP. The present discrepancy stands at 34(11) x 10 -10 with the above plots showing 
1 and 2 a contours. 
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Figure 3: Contours of the second generation sneutrino mass, m^, are plotted in the plane 
of (ttif, M\j2). The contours are in the units of GeV. The four plots, are obtained from 
the two minima, minimum A and minimum B with /i = 120 and 300 GeV as labelled. All 
points in the top left corner with approximately Mi/ 2 > 700 GeV and nip > 700 GeV are 
unphysical due to the lightest stau becoming the LSP. 
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Figure 4: This plot displays the effect on x 2 due to a future change in the value of the 
muon g — 2 discrepancy. The value of the muon anomalous magnetic moment is varied 
from the present value of 34 x 10 -10 down to zero. The resulting change in x 2 is observed 
for three points in the (mp, Mi/ 2) plane. 
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Figure 5: This plot displays the effect on x 2 due to an increase in the lower bound on the 
Higgs mass from direct searches. As in Figure 4 the variation in x 2 is observed through 
individual points in the (mp, Mi/ 2) plane. 
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Figure 6: Contours of sin(/3 — a), which defines the strength of the Z boson coupling to 
the Higgs h° relative to that of H°. For values of sin(/3 — a) near one the Z — h° coupling 
is large and for small values the Z — H° coupling is large. The contours are plotted using 
data from minimum A with /i = 120 GeV. The the best fit point at My^ = 450 GeV, 
mf = 500 GeV is marked with an asterisk. 
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Figure 7: Contours of BR(t — > ^7) are plotted in the plane of (thf, M]/ 2 ). The four 
plots, are obtained from the two minima, minimum A and minimum B with /i = 120 and 
300 GeV as labelled. All points in the top left corner with approximately M1/2 > 700 GeV 
and mp > 700 GeV are unphysical due to the lightest stau becoming the LSP. 
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Figure 8: Contours of r& = Yb/Y T are plotted in the plane of {mp, M1/2). The four 
plots, are obtained from the two minima, minimum A and minimum B with /i = 120 and 
300 GeV as labelled. All points in the top left corner with approximately M1/2 > 700 GeV 
and mp > 700 GeV are unphysical due to the lightest stau becoming the LSP. 



35 



M 1/2 [GeV] 



1000 

900 

800 

700 

600 

500 

400 
1000 

900 

800 

700 

600 

500 

400 



(i=120GeV 
" Min A 



0.95 



/ 1.07 
1.05 



1 .0 1.02 



1 1 1 — H- 

(i=120GeV 

Min B 



1.07 



1.05 
1.02 



1.0 



J I I I L 



(i=300GeV / 
' Min A 



0.9 



; , 1 .07 
1.05 



0.95 



1.0 



1.02 



(i=300GeV, 



Min B 



0.9 



1.07 



0.95 



, 1.05 
1.02 



1.0 

J I I I L 



400 500 600 700 800 9001000 500 600 700 800 9001000 

m F [GeV] 



Figure 9: Contours of rt = Yt/Y T are plotted in the plane of {rap^Mi^)- The four 
plots, are obtained from the two minima, minimum A and minimum B with fj, = 120 and 
300 GeV as labelled. All points in the top left corner with approximately Mi/ 2 > 700 GeV 
and mp > 700 GeV are unphysical due to the lightest stau becoming the LSP. 
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Figure 10: x 2 contours in the r t — plane. The plots are generated with fx = 120 GeV 
and for minimum A. The three plots are each generated with fixed M1/2, mp as labelled. 
The plots display the x 2 penalty which is required for exact Yukawa unification to be 
achieved. 
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Figure 11: Contours of the CP odd Pseudoscalar Higgs mass, m A o, are plotted in the 
plane of (mpjM^). The contours are in the units of GeV. The four plots, are obtained 
from the two minima, minimum A and minimum B with = 120 and 300 GeV as labelled. 
All points in the top left corner with approximately Mi/ 2 > 700 GeV and mp > 700 GeV 
are unphysical due to the lightest stau becoming the LSP. 
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Figure 12: The four panels contain scatter plots of the values of Am\ tm against sin 2 2#i3 
coming from the best fit points in the (tuf, Min) plane. Each plot shows results obtained 
from either minimum A or minimum B, with [i = 120 or 300 GeV as labelled. The points 
are grouped according to their \ 2 values as inticated. 
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